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Abstract 



In this work we construct algebraic equation for elements of spectrum of linearized Fokker - Planck 
differential operator for incompressible fluid. We calculate roots of this equation using simple numeric 
method. For all these roots real part is positive, that is corresponding solutions are damping. Eigen- 
functions of linearized Fokker - Planck differential operator for incompressible fluid are expressed as 
linear combinations of eigenfunctions of usual Fokker - Planck differential operator. Poisson's equation 
for pressure is derived from incompressibility condition. It is stated, that the pressure could be totally 
eliminated from dynamics equations. The Cauchy problem setup and solution method is presented. 
The role of zero pressure solutions as eigenfunctions for confluent eigenvalues is emphasized. 
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1 Introduction 

In our previous work [1] we derived linearized Fokker - Planck equation for incompressible fluid 
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where 



n = n{t,xi,X2,X3,vi,V2,vs) - density; 
p = p{t, xi, X2, a^s) - pressure; 
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t - time variable; 

xi,X2,xs - space coordinates; 

vi,V2,vz - velocities; 

a - coefficient of damping; 

k - coefficient of diffusion. 

No attempt was made to solve this equation. 

Peculiarity of this equation consists in the fact, that for two unknown variables n and p we have only 
one differential equation. This is enough, because there is additional normalization requirement (1) on n 
variable and p variable depends only on space coordinates x, y, z and time t. 

In [2] some simple solutions of nonlinear equation were studied. Particularly flow with zero pressure 
is of interest for our present studies, because for this flow the source of nonlinearity is absent. This flow 
solves both nonlinear and linear equations. 

Present work is devoted to solution of linearized equation. We consider only the case of parallelepiped 
with opposite sides identified (i.e. "periodic boundary conditions"). We try to formulate and solve 
Cauchy problem for linearized equation. 

2 Fourier decomposition of solution 

We know the form of Cauchy problem solution for the case of usual Fokker-Planck equation (see [3]). 

+OC +00 +00 oc oo oo 

(3) 

mi=— CO m2=— oo m3=— oo pi=0 P2=0p3=0 

where eigenfunctions of usual Fokker - Planck operator are 



cPm.m2msn.n2ns = n exp {2m^ix, - exp (-^^|) H^^ (^v, + ^^^j) • (4) 

It is only natural to seek solution of the present problem in the same form. We need only add 
expression for the new p variable 
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We try to represent the coefficient before in equation (2) in the same way as a sum of Fourier 



series 
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In expressions (3, 5, 6) we introduced following Fourier coefficients: 
{t) - coefficients of decomposition unknown variable n; 



mim2m3pip2P3 
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[t) - coefficients of decomposition of unknown variable p; 



(k) - coefficients of decomposition of known variable, which represent coefficient before p 



gradient. Values of B are presented below. 

Using these coefficients, we rewrite equation (2) as 
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Wc reduced the partial differential equation (2) to system of ordinary differential equations for Fourier 
coefficients. To proceed with solution, we need expressions for known coefficients Bmim2m3nin2n3{k) ■ 



3 Expressions for Fourier coefficients 



In this auxiliary section we find explicit expressions for Fourier coefficients of some known functions of 
velocities. These functions are products of multipliers, which depend only on one independent variable. 
Therefore we can consider only one velocity variable in this section. 

We start from definitions 
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Functions and V'mn are of course orthogonal (see [3]) 
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Let us find Fourier coefficients for following functions: 
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To find these coefficients, we need to calculate integrals 
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We calculated these integrals in our previous work [2] 
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Thus we get following expressions for Fourier coefficients 
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4 Dynamics of Fourier coefficients 

In this section we substitute expressions for known coefficients in terms of am2n2 and bmim coefficients, 
which we find in the last section, to the main equation (7). Namely, we use expressions 
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Then (7) reads 
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This is main equation, which describe dynamics of Fourier coefficients. We delay actual substitution 

of Q-mirij and bmiUj 

until (31). 



5 Incompressibility condition 



We use the incompressibility condition (1) to eliminate Pmim2m3 

(t) from (23). 
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So incompressibility condition is equivalent to following equation 
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Let us suppose, that coefficients A^^j^^^^nin2n3it) satisfy (26) at the moment t. They must satisfy 
this equation at the next moment t+dt. Values of A at the next moment are defined by dynamics equation 
(23), which contains besides A{t) also pressure P. Therefore incompressibility condition, written for the 
next moment, will give us equation for pressure P. Derivation of this equation is rather long procedure, 
which ends in equation (42). 
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Differentiate (26) on time t and get 
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We can easily calculate the sums on nj in RHS of (29). For this purpose let us introduce the 
sum on the group of terms with constant {ni + n2 + ns) = J. 



If) fe)" 



^ ^ Cmj^jC^27l2'^3773 ^ 

771+772+713 = '^ 



^On7]^jj]^077j2772^7773 773 



"T-l 6mi77i ^2 ^7(72772 ^3 ^m-. 



+ 



7^,3773 



Ol Omim O2 Of7l2772 ^3 O^gjjg 

Let us substitute to (30) values of coefficients Ojj , hij , Cij 
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The last sum according to Newton's binomial theorem is 
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We see, that result Sj depends on two variables J = ni + n2 + ns and M 
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With this variables (32) reads 
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This result is a little disappointment, because wc still not get desired equation for 
need insistence to achieve success. The result is already near. 
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Let us use (37) to remove terms with { ^^^^ ^ in (39). To calculate coefficient before -Pmim2m3(i) we 
perform once again summation on group of terms with constant {n\ + n2 + na) = J. The sum for each 
group is calculated as before, but this time this sum is multiplied by J before final summation on J 
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(42) is algebraic equation, but in the same time it is Fourier transform of some differential cqTiation 
for original unknown variable n. Namely according to (33) M = ((^^)^ + (^^)^ + (^^)^ 
. Each m? term in (42) is Fourier transform of second partial derivative of p on corresponding space 
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variable, their sum is Fourier transform of Laplace operator. Therefore (42) is Fourier transform of 
Poisson's equation for pressure. 

Solve this equation for Pmim2m3{t) and get 
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To get final result, substitute values of amn from (17) 
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We totally eliminated pressure from dynamic equation. Thus the problem is reduced to the system of 
ordinary linear differential equations (46). Such systems are solved by Euler's exponential substitution. 
The only difficulty is, that the number of variables is indefinitely great. We try to integrate (46) using 
special properties of its matrix. 
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6 Special eigenvalue problem 

In this section we consider purely algebraic eigenvalue problem for special form of matrices. 
Let the matrix have the following form 

Mij = Dij + P,j; (47) 

where Dij - diagonal matrix 

D^, = {': (48) 

ydi- 1 = J 

Pij - diadic product matrix 

Pij = krj. (49) 

The eigenvalue problem consists of following tasks: 

1) To find a set of eigenvectors x such, that matrix product of M and x is proportional to x 

MijXi = Xxj. (50) 

2) To find corresponding set of proportionality coefficients - eigenvalues A. 

For our special form of matrix 
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This gives very simple expression for components of eigenvector, provided eigenvalue A is evident 
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Substitute this expression to (52) and get 



Two cases are possible here. The common case is S\ = 0. 

(56) gives equation for A. When all di are different, algebraic equation (56) has degree n, where n - 
dimension of matrix M. 



10 



If all roots of (56) are different, we get the full set of eigenvalues, then from (54) we find full set of 
eigenvectors. Exact value of S is of no meaning, we can put for example S = 1 in (54). 

For example, when all {Ikrk) are positive or all are negative and all di are real (and different - see 
above), one could guarantee that all n roots of (56) are real and different. This follows from the fact, that 
di separate the roots of (56). Therefore there are n — 1 roots between di and one more root ximin{di) 
(the case kriO). This fact makes numeric evaluation of roots rather simple. Unfortunately we deal with 
quite opposite case - signs of our {Ikrk) are alternating. Nevertheless we shall find roots - see below. 

When some di are equal, equation (56) has less roots then matrix degree. We can consider this case as 
confluent. In this case additional eigenvalues A (besides roots of (56)) must be equal to iterated value di 
and we must put for these eigenvalues S* = in (53). Then all Xi , besides that in columns, corresponding 
to iterated di, must be zero. The rest nonzero Xi must satisfy orthogonality condition (52) (with S = 0). 

Let us consider the eigenvalue problem for transposed matrix . Let us denote y - eigenvectors for 
M^. Then 

{Dij - A dij)yi + {lkyk)rj = 0; (57) 
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We get for A equation (56) once again - eigenvalues of conjugated problems are equal. Components 
of eigenvectors for conjugated problem are calculated from (59), where S is arbitrary nonzero number, 
for example S = 1. 

Eigenvector for conjugated problems with different eigenvalues A and fx are orthogonal: 
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Let us return to our equation (46). Comparing with previous section, we could identify matrix components 
in the following way. 
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case - see previous section. Some eigenvalues are equal to c?nin2Ti3 • 
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Nondiagonal components of the matrix are equal to: 
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Let us once again calculate partial sum Sj of terms, for which ni + r^2 + 713 = J. We have 
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We see, that all eigenvalues A are proportional to 
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where coefficients ^ are roots of equation 
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M + ^— M-(^-«) [7((2 + M - 0, M) - 7((1 + M - 0, M)] = 0. (71) 



We can further simplify (71) using known recurrence (see [4], 9.2) 
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One exact root of equation (75) we can find easily - this is root 
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Another form of (77) we obtain using modified incomplete gamma function 
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1*iP,x) = Y^j{p,x). (78) 

This form is 

1 + Ce^r(l + M - 7*((1 + M - 0,M) = 0. (79) 

Advantage of this form is that 'y*{p,x) is a single valued analytic function of p and x possessing no 
finite singularities. 

We calculate some roots of equation (77) using definition of incomplete gamma-function (66 - 67). 
Namely we keep only finite number of terms in the sum (66 - 67) and solve resulting algebraic equation 
using Newton's method. As the number of terms increase, the roots converge rapidly. 

As initial approximation for roots we use position of poles, that is we search roots in the close vicinity 
of the pole. We keep corresponding term and approximate contribution of the rest poles by two first 
terms of Taylor series. So we have quadratic equation. When this equation has two conjugated roots, 
Newton's iterations converge to two conjugated roots of full equation. When quadratic equation has two 
real roots, one root is really located in the poles vicinity and another is far enough, so that Newton's 
iterations diverge. 
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The results of our calculations are presented in the Table 1 (see APPENDIX). We see, that: 

- All roots have positive real part. Therefore exists root with the least real part and roots can be ordered 
according their real part in ascending order. In the following we suppose that such ordering is done. 

- For each M only finite number of roots are complex and imaginary part of roots decreases with root 
number. 

- Starting from some root all roots are real (tail). 

- Real roots are located very close to poles (natural numbers) and rapidly get indistinguishable. 

On this stage we shall content ourself by these experimental results and shall not try to give them 
rigorous proof. The theory of distribution of roots of analytic functions in question is rather ample (see 
[6], [7], [8]). 

As we see from previous section, besides roots of (77) there exist eigenvalues, which are exactly 
equal to diagonal values (inin2n3 (62). They correspond to zero-pressure solutions from our work [2]. 
Orthogonality conditions (52), which must be satisfied by eigenfunctions, according to (64) and (43) 
mean, that Pmimzma (t) = 0. 



8 Solution of Cauchy problem 

In this section we describe construction of solution of linearized Fokker - Planck equation for incompress- 
ible fluid. 

1) To setup Cauchy problem we must set initial value of n variable. This initial value no = 
no{xi,X2,X2,vi,V2,V3) must satisfy incompressibility condition (1). There is no need to set initial value 
of pressure, because it is fully determined by n (see (42)). 

2) Calculate Fourier coefficients ^mim2m3pip2P3(0) (see [3] for details): 

/ a \ 2 1 



2Pi+P2+P3pi\p2]p^] \2TTk J 010203 



X exp 



2k V 

a'i foo rco roo 



/ dxi / dX2 dX3 / / no{xi,X2,X3,Vi,V2,V3) 1pmim2m3PiP2P3dvidV2dV3. 

Jo Jo Jo J -00 J -00 J -00 



where 

7/;^ = I I p-^n I — 9.Tri- , ^^n- w ' 

a a J ^ \ \ 2k \ a^a 



'(Pmim2m3nin2n3 = H ( -2vri^(a;j + -^) ) H^^ ( \/ ( + ^2^.^ ] ] ' (^1) 



3) Change from eigenfunctions of simple Fokker - Planck equations to eigenfunctions of linearized 
Fokker - Planck equation for incompressible fluid. According to (61) projection of vector on eigenvector 
x^, corresponding to eigenvalue A, is 
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or according to (64) 



*,I/lfc'2fc'3 



X 



/27rmiyi /27rm2y2 /2TTm?,y ^ 
\ aai ) V OLa2 J \ aa^ J 

(1/1 + 1/2 + ^3)^^ I '^nirt2n3 



(A dii^iir^vz) J \ it ^"^ di^ii^^i^s)^ J 

where dnmzns is defined by (62), r„j„2„3 is defined by (64). 

4) Initial field can contain some zero pressure solutions. Let us suppose, that there exist a group of 
coefficients with constant J = i/; , for which 

This group docs not contribute to any Amim2m3X (sec 83). Such groups, if present, we must consider 
separately. They correspond to zero pressure solutions of our work [2]. 

5) Given initial values of Ax we can calculate their values for the arbitrary moment t according to 
exponential law 

■Amim2m3x(j') ~ C ^mim2m3A(0)- (^^) 

6) Evolution of zero pressure solution is determined by exponential multiplier e~'*('^+^)*. 

7) Inverse transition from Ax to Aj. is 



or according to (63) 



A A ^^, exp(M) 

A 

1 /2l--\('^i+'^'^+'^s)/'^ 



(87) 



2"i+"-2+n3„^!„2!n3! V a. 



/27rmi\"i /27rm2\"2 /27rm3y3<^ 
\ aa J \ aa J \ aa J 



( lirnij 
arin + k ' 



X 

\ 2' 



(A dfi^ji2n3^ 



c^nmans defined by (62). 

8) Contribution from zero pressure solutions must be added to (87). 

9) Pressure for each moment of time is determined by (44) . 
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DISCUSSION 



We see that spectral properties of Fokker - Planck linearized differential operator for incompressible fluid 
are different from properties of usual operator. General spectrum structure is roughly similar, but nearest 
to zero eigenvalues are complex. Therefore most slowly damping modes are most strongly vibrating - 
very interesting result. Generally all modes are damping with time, flows tend to the rest. 
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APPENDIX 1 

Roots of equation (71) 



Root ^ = M is omitted 



NN 


M=l 


M=2 


M=3 


1 


3.84958810 - 1.92315575 i 


4.52745332 - 3.27206660 i 


5.04119504 - 4.34568739 i 


2 


3.84958810 + 1.92315575 i 


4.52745332 + 3.27206660 i 


5.04119504 + 4.34568739 i 


3 


5.94063198 - 1.14455587 i 


7.06814092 - 2.67480339 i 


7.90915778 - 3.92341797 i 


4 


5.94063198 + 1.14455587 i 


7.06814092 + 2.67480339 i 


7.90915778 + 3.92341797 i 


5 


7.69165960 - 0.29172595 i 


9.14966499 - 1.97226724 i 


10.24525406 - 3.33703003 i 


(3 


7.G91659G0 + 0.29172595 i 


9.149GG499 + 1.9722G721 i 


10.2-1525 lOG + 3.33703003 i 


7 


11.000932211 


11.00141191 - 1.21959567 i 


12.31747194 - 2.67114557 i 


8 


11.999889388 


11.00141191 + 1.21959567 i 


12.31747194 + 2.67114557 i 


9 


13.000011755 


12.72299911 - 0.42498830 i 


14.22286668 - 1.95786294 i 


10 


13.999998866 


12.72299911 + 0.42498830 i 


14.22286668 + 1.95786294 i 


11 


15.000000099 


14.07454792 


16.00982241 - 1.21325010 i 


12 


15.999999999 


14.98376860 


16.00982241 + 1.21325010 i 


13 


17.000000000 


16.00272440 


17.72151356 - 0.43464100 i 


14 


17.999999999 


16.99956146 


17.72151356 + 0.43464100 i 


15 


19.000000000 


18.00006506 


19.98044555 


16 


19.999999999 


18.99999098 


19.08192627 


17 


21.000000000 


20.00000116 


19.98044554 


18 


22.000000000 


20.99999985 


21.00363177 


19 


23.000000000 


22.00000001 


21.99933574 


20 


24.000000000 


22.99999999 


23.00011373 
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NN 


M=4 


M=5 


M=6 


1 


5.-17198-106 - 5.26768710 i 


5.85889775 - G. 08913226 i 


G. 20771808 - G. 83731700 i 


2 


5.47498406 + 5.26768710 i 


5.85889775 + 6.08913226 i 


6.20771808 + 6.83731700 i 


3 


8.61344667 - 5.01144505 i 


9.23406886 - 5.99054148 i 


9.79675137 - 6.88888105 i 


4 


8.61344667 + 5.01144505 i 


9.23406886 + 5.99054148 i 


9.79675137 + 6.88888105 i 


5 


11.15620212 - 4.53812831 i 


11.95509246 - 5.62665552 i 


12.67703265 - 6.63086922 i 


6 


11.15620212 + 4.53812831 i 


11.95509246 + 5.62665552 i 


12.67703265 + 6.63086922 i 


7 


13.40577133 - 3.95833350 i 


14.35634208 - 5.13111612 i 


15.21274497 - 6.21756080 i 


8 


13.40577133 + 3.95833350 i 


14.35634208 + 5.13111612 i 


15.21274497 + 6.21756080 i 


9 


15.47089306 - 3.31449740 i 


16.55737228 - 4.55585073 i 


17.53366817 - 5.70963951 i 


10 


15.47089306 + 3.31449740 i 


16.55737228 + 4.55585073 i 


17.53366817 + 5.70963951 i 


11 


17.40539050 - 2.62778318 i 


18.61704841 - 3.92690815 i 


19.70341924 - 5.13774003 i 


12 


17.40539050 + 2.62778318 i 


18.61704841 + 3.92690815 i 


19.70341924 + 5.13774003 i 


13 


19.24064677 - 1.91036878 i 


20.56955610 - 3.25932488 i 


21.75880445 - 4.51960861 i 


14 


19.24064677 + 1.91036878 i 


20.56955610 + 3.25932488 i 


21.75880445 + 4.51960861 i 


15 


20.99699656 - 1.17010968 i 


22.43688595 - 2.56260936 i 


23.72343062 - 3.86649261 i 


16 


20.99699656 + 1.17010968 i 


22.43688595 + 2.56260936 i 


23.72343062 + 3.86649261 i 


17 


22.70576400 - 0.39532166 i 


24.23420403 - 1.84317596 i 


25.61354419 - 3.18599046 i 


18 


22.70576400 + 0.39532166 i 


24.23420403 + 1.84317596 i 


25.61354419 + 3.18599046 i 


19 


24.07528691 


25.97268205 - 1.10579718 i 


27.44091387 - 2.48348667 i 


20 


24.98146932 


25.97268205 + 1.10579718 i 


27.44091387 + 2.48348667 i 



# 


M=7 


M=8 


M=9 


1 


6.53001156 - 7.52895368 i 


6.83127967 - 8.17518985 i 


7.11531330 - 8.78391538 i 


2 


6.53001156 + 7.52895368 i 


6.83127967 + 8.17518985 i 


7.11531330 + 8.78391538 i 


3 


10.31618762 - 7.72399088 i 


10.80167947 - 8.50772212 i 


11.25955129 - 9.24856515 i 


4 


10.31618762 + 7.72399088 i 


10.80167947 + 8.50772212 i 


11.25955129 + 9.24856515 i 


5 


13.34199931 - 7.56851364 i 


13.96257983 - 8.45167276 i 


14.54728640 - 9.28905550 i 


6 


13.34199931 + 7.56851364 i 


13.96257983 + 8.45167276 i 


14.54728640 + 9.28905550 i 


7 


15.99974275 - 7.23543754 i 


16.73291571 - 8.19691791 i 


17.42277980 - 9.11081249 i 


8 


15.99974275 + 7.23543754 i 


16.73291571 + 8.19691791 i 


17.42277980 + 9.11081249 i 


9 


18.42897216 - 6.79355442 i 


19.26164201 - 7.81977962 i 


20.04405940 - 8.79717049 i 


10 


18.42897216 + 6.79355442 i 


19.26164201 + 7.81977962 i 


20.04405940 + 8.79717049 i 


11 


20.69785002 - 6.27781735 i 


21.62130874 - 7.35928644 i 


22.48793355 - 8.39101050 i 


12 


20.69785002 + 6.27781735 i 


21.62130874 + 7.35928644 i 


22.48793355 + 8.39101050 i 


13 


22.84567896 - 5.70854973 i 


23.85362358 - 6.83821861 i 


24.79844419 - 7.91745869 i 


14 


22.84567896 + 5.70854973 i 


23.85362358 + 6.83821861 i 


24.79844419 + 7.91745869 i 


15 


24.89759650 - 5.09866262 i 


25.98519281 - 6.27108315 i 


27.00361964 - 7.39255586 i 


16 


24.89759650 + 5.09866262 i 


25.98519281 + 6.27108315 i 


27.00361964 + 7.39255586 i 


17 


26.87086392 - 4.45689278 i 


28.03424697 - 5.66771323 i 


29.12262537 - 6.82719343 i 


18 


26.87086392 + 4.45689278 i 


28.03424697 + 5.66771323 i 


29.12262537 + 6.82719343 i 


19 


28.77796298 - 3.78943945 i 


30.01394562 - 5.03509475 i 


31.16927222 - 6.22911873 i 


20 


28.77796298 + 3.78943945 i 


30.01394562 + 5.03509475 i 


31.16927222 + 6.22911873 i 
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